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We have begun a series of experiments on mixed bosonic quantum fluids. 
Our system is mixed Bose-Einstein condensates in dilute Rb-81. By simul- 
taneously trapping the atoms in two different hyperfine states, we are able 
to study the dynamics of component separation and of the relative quantum 
phase of two interpenetrating condensates. Population can be converted from 
one state to the other at a rate that is sensitive to the relative quantum phase. 

PACS numbers: 03.75. Fi, 05.30.Jp, 32.80.Pj, 42.50. Dv 



1. INTRODUCTION 



The observation of Bose-Einstein condensation (BEC) in a dilute gasl 




has led to a new generation of experiments in quantum fluids. Historically, 
in the study of quantum fluids, some of the most intriguing behavior has 
been found in the behavior of fluid mixtures. We have thus been motivated 
to pursue a series of mixed-condensate experiments in Rubidium-87.@~@ The 
experimental tools and the theoretical world view of atomic physics are al- 
most entirely disjoint from those of traditional low-temperature physics, and 
this lack of overlap has led to some confusion. The purpose of this paper 
is to review some of our early results on mixed condensates, and to provide 
a qualitative exegesis of the theoretical and experimental techniques that 
underpin our work. Of necessity we have omitted here much of the detail 
provided in the original papers. oQ The MIT BEC group also has begun 
experiments in mixed condensates B 

The experiments described here all begin with a sample of approxi- 
mately 5 x 10 5 spin-aligned, Bose-condensed Rb-87 atoms at a temperature 
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of less than 50 nK, confined in a 3-d, harmonic magnetic potential at density 
1 x 10 /cm 3 . The number density of noncondensed atoms is a factor of 30 
or more smaller. In this paper we will skip the details of refrigeration — the 
technique is a hybrid method that combines many of the notable advances^ 
in atom cooling from the 1980s: laser cooling and trapping, magnetic trap- 
ping, and evaporative cooling. The synthesis of these various approaches 
was begun at JILA from 1989-1992,BEl but it was only after considerable 
additional developmentE^Ej that, in 1995, the parts worked well enough to- 
gether to permit the first observation of BEC in dilute atomic gases.Era In 
the current version of our apparatus, the cooling process is automated and 
reliable, much like a (well-behaved) dilution refrigerator. A freshly prepared 
condensate is available for study about once per minute. The isolation of 
the sample from the 300 K environment (less than 1 cm away) is excellent 
- thermal equilibration times with the environment are on the order of 
centuries — but residual heating and the collisional decay of the conden- 
sate limit the duration of an experiment on any given condensate to a few 
seconds. 

The ground state of an Rb-87 atom is split in a magnetic field into eight 
levels by hyperfine and Zeeman interactions (Fig. |]). The states labelled |1), 
1 2), and |3) can be trapped at a local minimum in the magnetic field; it is 
with states |1) and |2) (\F = l,m = —1), and \F = 2,m = IV. respectively) 
that we have performed most of our mixed-fluids experiments.O'E3 The other 
relevant internal structure of the Rb-87 atom is a set of excited-state levels 
1.6 eV above the ground state (Fig. [l]). We probe the spatial distribution of 
the atoms by imaging the absorption on an optical transition to these levels. 
The natural 6 MHz linewidth of the optical transition is much less than the 
energy splitting between, for instance, states |1) and |2). With the correct 
choice of laser frequency, we can image the spatial distribution of either state 
|1) or of state |2), or, if we prefer, of the combined density of the two states.! 

The paradoxical aspect of the Rb-87 mixed-condensate system is that 
it combines two traits that are usually thought of as mutually exclusive: 
(i) the |1) and |2) components of the condensate are distinguishable, in the 
usual senses of the term — the components do not at experimental energy 
scales spontaneously interconvert, and the components can be selectively 
detected and imaged; and (ii) the components possess a measurable relative 
quantum phase, which evolves at a rate proportional to the difference in 
chemical potential between the two fluids, and which can, under the right 
circumstances, be read out, much as the current across a Josephson junction 
reads out the relative quantum phase across it. Interconversion (and phase 
read-out) can be driven via a stimulated process and therefore can be turned 
on and off. The system thus bears considerable resemblance to the idealized 
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5P 3/2 excited states 




F = 2 



m = -2 m = -1 m = m = +1 m = +2 

Fig. 1. The relevant hyperfine and Zeeman structure of Rb-87 in a magnetic 
field. The atomsare moved coherently from state |1) to state |2) with a 
microwave pulse. £9 Imaging is accomplished via absorption in optical tran- 
sitions to the excited states. 

pair of condensates linked by a removable Josephson junction that LeggettS 
and others have envisioned in various thought experiments. 

In Section 2 we will discuss distinguishable-fluid experiments; in Sec- 
tion 3 we discuss simple phase-coherent experiments; in Section 4 we discuss 
experiments that combine the two aspects; and in Section 5 we mention 
various future experiments with weak coupling. 



2. INTERPENETRATING, DISTINGUISHABLE FLUIDS 

Species |1) and |2) meet both usual criteria for calling fluids "distin- 
guishable." One can in the literal sense distinguish one density distribution 
from the other, simply by adjusting the frequency of the probe laser so 
that one species or the other casts its shadow onto our imaging CCD array 
(Fig. |). This measurement, as performed in Refs. ||-[7|, is destructive, but as 
has been shownS it needn't be so. In any case the fact that we can prepare 
condensates with a high repetition rate and a large degree of reproducibil- 
ity means that the destructive aspect of the imaging is not very important 
to us: we observe the time-evolution of the clouds simply by repeating the 
experiment many times with increasing dwell-times. The second sense of 
"distinguishable" is that particles do not spontaneously interconvert. This 
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requirement is enforced in our system by the enormous difference in internal 
energies between the two states. In the absence of applied electromagnetic 
fields, the 6.8 GHz hyperfine energy (Fig. [I]) is a million times larger than 
any other energy scale of the condensate system. The energy released by 
a single atom converting from state |2) to state |1), if that energy were to 
be distributed thermally through our sample, is sufficient to drive the entire 
sample out of the Bose-condensed state. The fact that we don't observe our 
condensates melting during our measurements is a guarantee that sponta- 
neous inter conversion is not a factor in the system. 

The positions of the two fluids are determined by the their respective 
confining magnetic potentials. The magnetic moments of the two states are 
nominally the same, but due to various small effects — gravity, the nuclear 
magnetic moment the onset of nonlinearity in the Zeeman shifts, and subtle 
dynamical effects"'" of our magnetic trap — the location of the two min- 
ima of the confining potentials V\(r) and l^(r) can be adjusted to be either 
exactly coincident or slightly offset. For the work described here, the con- 
fining potentials are axially symmetric and harmonic, with single-particle 
oscillation frequencies of uj z /2tt = 60 Hz, and u) t /2 , k = 21 Hz. With the 
spatial offset between the two potentials set to zero, and in the absence of 
interspecies interactions, the two fluids would be completely interpenetrat- 
ing. In fact, inter- (and intra-) species interactions are a dominant factor in 
determining the density distributions. The dynamics are well-described in 
the mean-field language of coupled Gross-Pitaevskii equations for the order 
parameters:B3 

ih®x = ( ~^V 2 + Vi + uil^l 2 + ni 2 |$ 2 | 2 J $i (1) 



and 



ih® 2 = ( ~ J^ y2 + ^2 + Vm + u 2 \®2\ 2 + ^2i|^i| 2 j $2 (2) 



where m is the mass of the Rb atom, Vm is the magnetic field-dependent 
hyperfine splitting between the two states in the absence of interactions, 
Ui = 4:7rh 2 ai/m and mj = Airh 2 aij/m, |$j| 2 is the condensate density, and 
the intraspecies and interspecies scattering lengths are a, and a^- = Ojj. Note 
that these equations do not allow for species inter conversion; they conserve 
independently the number of atoms in each species. 



Consider, 
state solutions^ 2 . 



e theoretical work has already gone into studying ground- 
and small- amplitude excitationseJ of the order parameters. 
We review here the results most relevant to the current experiments. As 
the total number of atoms trapped increases, the quantum kinetic energy 
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term V 2 $ becomes less significant, until it is almost negligible for the num- 
bers used in double-condensate experiments. Most theoretical treatments 
neglect this term, thus making the so-called Thomas-Fermi approximation. 
In this case, the steady-state, single-component condensate density profile 
in a parabolic potential will have an inverted parabaloid density profile, 
with density tapering smoothly from theneak at cloud center to zero at its 
edge. Due to an accidental degeneracyBo in the Rb-87 scattering system, 
a\ as ayi « a,2- For this special case, one can see by inspection of Eqs. || and § 
that the steady-state total density of a two-component BEC will also have 
an inverted parabola form, even though the relative densities of the compo- 
nents may have intricate structure. It turns out there is a critical valuel! for 
the interaction term a\ 2 = y^l^i- For ai 2 > af 2 , two components should 
have little spatial overlap — they should spontaneously separate in the trap, 
while for a 12 < af 2 , there should be a relatively large region of interpenetra- 
tion. In Rb-87, the scattering lengths are known at the 1% level to be in the 
proportion ai : a\2 : 02 - 1-03 : 1 : 0.97, with the average of the three being 
55(3) A.IHi 

Since the mutual interaction term is within 1% of its critical value, any- 
thing that breaks the symmetry between the two species becomes important. 
For instance, the fact that a± > a 2 means it is energetically favorable for the 
|1) atoms to preferentially move towards the lower density region at the pe- 
riphery of the cloud, forming a spherical shell around |2),E3 but since the 
difference between a\ and a 2 is small, even a rather minor vertical offset in 
the spatial centers of V\ and V2 will result in the components' separating 
up-and-down, rather than radially in-and-out. In steady-state, then, one ex- 
pects the clouds to be largely, although not completelynj spatially separated. 
Further, because the steady-state energetics only modestly favor component 
separation, one would expect that in the ensuingdynamical behavior, small 
oscillations about the steady-state configurations will feature a number of 
"soft" modes, with frequencies low compared to the excitations of the total 
density. 

All of these qualitative theoretical expectations have been borne out 
in preliminary experiments.! Figure || shows the density profiled of two 
condensates after they have come to steady-state. In this measurement, the 
potentials have been offset a distance equal to only 3% of the overall sample 
size, but the resulting component separation is much larger. Note however 
that there does remain a region of overlap, which is useful for work described 
below. The |2) component is originally created (as described in Section 3.2, 
below) such that it completely overlaps the |1) component, and we have 
observed! the time evolution of the component separation as it relaxes to 
the steady-state, separated condition. During this time, the overall density 
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Fig. 2. Contour plot of the expandedlll density distribution after the two 
components have undergone separation and their relative motions have 
damped. The image is 182 microns square. Before expansion the size of 
the total distribution was 40 by 14 microns. Note that a considerable region 
of interpenetration remains. 



profile is essentially unperturbed. The period of the resulting damped oscil- 
lations of the relative positions of the two components is about 30 ms. In 
comparison, the period of the lowest-order nontrivial excitation of the total 
density profile is about 25 ms. 

If we do not deliberately break symmetry of the confining potentials in 
the experiment but rather keep them as perfectly overlapped as possible, 
we see, at early times, the onset of the expected inward-outward component 
separation. At longer times the inner ball of |2) atoms does not stay well- 
centered in the hollow shell of |1) atoms, but drifts horizontally towards a 
preferred side. There appears to be an as-yet uncharacterized imperfection in 
the trap that breaks the radial symmetry of the experimental environment. 

The Rubidium-87 mixture is well-suited to a variety of experiments in 
mixtures of dilute bose gases. There are a host of experimental parameters 
to explore: the radiakaxial aspect ratio of the confining potential (from 3:1 
to greater than 1:30, depending on the form of magnetic trap used); the 
relative number in the two components; and the temperature. The critical 
temperature of a particular component depends on the number of atoms in 
that state. Thus, one can arrange for one component to remain in the con- 
densed state while the other one converts from the normal to the condensed 
state or vice-versa. All of these changes will have profound effects on the 
steady-state configuration and the spectrum of collective excitations. For 
example, imagine an experiment to study the effects of finite temperature 
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on the excitations.Ea At finite temperature, there will be in effect four differ- 
ent co-trapped fluids: normal and condensate fractions of both components. 
The normal fraction of a given component will find the condensate fraction 
of its own component twice as repulsive as it finds the condensate fraction 
of a dissimilar component (due to exchange terms.) Thus if one chooses a 
potential in which condensate |1) and condensate |2) tend to spatially sep- 
arate, one can have a situation in which the normal fraction of component 
|1) tends to preferentially collect at the spatial location of condensate |2), 
and vice-versa. What will the collective excitation spectrum look like? If 
we were to describe these excitations with the word "sound," what ordinal 
number should modify it? 



3. PHASE COHERENCE 

In the absence of any applied coupling the 6.8 GHz energy separation be- 
tween the two states is enormous, but the presence of a oscillating magnetic 
field tuned close to the hyperfine splitting changes the situation dramati- 
cally. An atom can readily be transferred from one state to the other - 
the difference in energy is absorbed from (or stimulated into) the coherent 
field of 6.8 GHz photons.c3 That the population transfer resulting from the 
application of a coupling field is sensitive to the relative phase of the two 
states is well-known in microwave spectroscopy^ and NMR and is easy to 
understand in the case of a single, two-level atom. 



3.1. Single-Atom Case 

The Schrodinger equation^ for the internal state of a two-level atom 
(or of any two-level quantum system) in the presence of a coupling drive at 
frequency u T f is 

%A X = cMi + ^e^A 2 (3) 

iA 2 = u 2 A 2 + ^-e-^Ai (4) 

where Ai is the probability amplitudes to be in state \i), ftwi is the internal 
energy of state \i), and f2 is the amplitude of the coupling. We assume the 
detuning 5 = u> r f — uj 2 + u>i is small, such that \5\ <C w r f , and concentrate on 
the case that the coupling drive is off (Q = 0) most of the time and on only 
for brief pulses of duration r, such that 1/r 3> \5\. When we discuss "phase" 
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in this context, we are talking about the real quantities a±, a 2 , or a r f defined 
in the relations A x = \A x \e iai , A 2 = \A 2 \e ici2 , and Ue'^ = |Sl|e ia rf. 

There is one particularly useful special-case solution of Eqs. || and f|. A 
pulse with duration and amplitude such that |f2|r = ^ is known as a 
pulse." If the atom is initially in state |1) with unit probability (i.e., |^4i| = 1 
and \A 2 \ = 0), then after the application of the pulse the atom will have equal 
probability of being in either state (i.e., \A\\ = l/\/2and \A 2 \ = l/\/2)- The 
^-pulse thus creates a 50-50 coherent superposition of the two states, also 
"writing" a particular initial relative phase. After the pulse is applied, the 
relative phase evolves at the rate uj 2 — u>\. The accumulated phase can be 
measured by a second ^-pulse.Ej When a ^-pulse is applied to a two-level 
system in 50-50 coherent superposition, the population transferred by the 
pulse from state 1 1) to state |2) is proportional to cos(ai — a 2 — a r f)> where 
the phases a refer to the total phase accumulated during the dwell time T 
between the two pulses. One measures the transition probability from state 1 
to state 2 by repeating the experiment for different T; the result is "Ramsey 
fringes," a cosine at the detuning frequency 5. This technique is known as 
the method of separated oscillatory fieldsE3 or "twin-pulse" spectroscopy. 



3.2. Application to Order Parameters 

Now it becomes clear how we modify Eqs. [I] and ^ to include a coupling 
field applied to a bose condensate: 



ih^ = (-^V 2 + 14+n 1 |^ 1 | 2 + u 12 \$ 2 \ 2 ) *i + ^Pe^** 2 (5) 
\ 2m I 2 

and 

im 2 = ( -|-v 2 + v 2 + v h{ + -u 2 |$ 2 | 2 + u 21 |$!| 2 ) $ 2 + ^l e -^ f *$ 1 

\ 2m J 2 

(6) 

The probability amplitudes Ai of Eqs. ^ and |I| have now become field ampli- 
tudes <&i, and the physics has become correspondingly richer. One can still 
perform experiments, however, that explore simple limits. 

In one such experiment, we begin with atoms in the pure |1), relaxed 
to a steady-state, near-pure condensate with spatial distribution described 
by $o- Next we apply a Impulse, transferring half the population to state 
1 2). Immediately after the |-pulse, the two-state system can be described 
by $1 = $ e iMl VV2> and $2 = Qoe^/y/ 7 }, where m and /i 2 , which differ 
by about the hyperfine frequency, are the chemical potentials of the two 
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Fig. 3. Population transfer from |1) to \2) resulting from twin Impulse cou- 
pling pulses, as a function of delay between the two pulses. The first pulse 
prepares the condensate in an equal superposition of the two states, the sec- 
ond pulse induces further population transfer sensitive to the relative phase 
that has evolved during the pulse separation time. The coupling drive is 
detuned from the energy difference between the two condensates by about 
360 Hz. 



components. We know however that this perfectly overlapping state is not 
the steady-state of the mixed BEC system; this is in fact the point of depar- 
ture for the relaxation-to-equilibrium experiments described in Section 2.1. 
Within a matter of a few tens of milliseconds the components will separate, 
but on shorter time scales, the condensates don't have time to realize that 
they are a pair of mutually repulsive quantum fluids. One could equally well 
describe the system at short times as a single condensate containing about 
a million atoms, each in an internal, coherent superposition governed by 
Eqs. | and |. Fi gure y shows the result of applying a second ^-pulse imme- 
diately, and then measuring the final number of atoms in the \2) state. The 
high-contrast Ramsey fringes that result are much as one would expect from 
a similar experiment performed on, for instance, an atomic beam.@ This 
"conventional" response of condensates at short times to coupling drives has 
also been seen in condensate Rabi oscillations 



4. COHERENT FLUIDS 

If the system is allowed to evolve longer after the first ^-pulse, so that 
the components begin to separate, it is no longer very useful to describe the 
system as a single condensate of atoms in a superposition state. One should 
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Fig. 4. Population transfer from |1) to \2) resulting from twin Impulse cou- 
pling pulses, as a function of delay between the two pulses. Results are 
similar to those shown in Fig. ||], except that over the longer duration the 
components physically separate — the loss of spatial overlap reduces the 
contrast ratio of the Ramsey fringes. 

think of them instead as distinct fluids evolving their separate ways, albeit 
(as we shall see) with a well-defined relative phase. The effects of compo- 
nent separation show up clearly in twin-pulse spectroscopy: Figure ^| shows 
some condensate Ramsey fringes collected over longer periods. The loss of 
contrast in the fringes is not due to inhomogenous broadening of the atomic 
transition, at least not in the conventional sense. Rather, it arises from the 
dwindling spatial overlap between the two components, which in turn arises 
from the mutual repulsion of the components, as we saw in Section 2. Dur- 
ing component separation the components are of course moving with respect 
to each other, and therefore there is a gradient in the relative phase across 
the sample; it is no longer so simple to define a single quantity as the phase 
differ erence between the two clouds. 

At still longer inter-pulse times, the results of the twin-pulse experi- 
ments become once again explainable in terms of a single relative phase. 
After the component-separation process has gone to completion, the relative 
motion of the condensates damps, and presumably the gradients of the rel- 
ative phase vanish as well. As shown in Fig. ^ above, the post-separation 
clouds preserve a reasonable amount of spatial overlap. As shown schemat- 
ically in Fig. ||, the application the second ^-pulse pulse at thi& time acts 
as sort of recombiner on a separated-arm atom-interferometer.Ea Depend- 
ing on the relative phase between the two components, we can see (Fig. |6|) 
destructive or constructive interference in the region where the components 
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expansion 




Fig. 5. A schematic^ of the condensate interferometer, (a) The experiment 
begins with all of the atoms in condensate |1) at steady-state, (b) After 
the first ^-pulse, the condensate has been split into two components with a 
well-defined initial relative phase, (c) The components begin to separate in a 
complicated fashion due to mutual repulsion as well as a 0.4 /jm vertical offset 
in the confining potentials (see also Fig. 3 of Ref. ||). (d) The relative motion 
between the components eventually damps with the clouds mutually offset 
but with some residual overlap. Relative phase continues to accumulate 
between the condensates until (e) at time T a second ^-pulse remixes the 
components; the two possible paths by which the condensate can arrive in 
one of the two states in the hatched regions interfere, (f ) The cloud is released 
immediately after the second pulse and allowed to expand for imaging. In 
the case shown, the relative phase between the two states at the time of the 
second pulse was such as to lead to destructive interference in the 1 1} state 
and a corresponding constructive interference in the |2) state. 
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Fig. 6. The value of the condensate density in the |2) state is extracted 
at the center of the overlap region (inset) and plotted as a function of T. 
Each point represents the average of 4 separate realizations and the thin 
bars denote the rms scatter in the measured interference for an individual 
realization. The thick lines are sinusoidal fits to the data, from which we 
extract the angular frequency \ii — [i\ — uj T f . 

overlap. Thus, the double-pulse experiment with a long-time between the 
pulses acts as a condensate interferometer that allows us to determine the 
relative phase between two largely separated condensate components. 

We find the contrast ratio manifest in Fig. |6| quite remarkable. It is 
evidence that the two components have a well-preserved memory of their 
relative phase even after there has been damping of the external degrees of 
freedom, with which the internal states are entangled. The phase of the con- 
densate evidently has a robustness one would not expect in a single-particle 
experiment. This is clearly an area worthy of extensive study - how robust 
is the relative phase, as a function of temperature, of relative equilibrium 
separation, of time between pulses? How well can the relative phase survive 
a "nondestructive" observation ,E9 which can take many sequential images of 
the same condensate? In single-atom interferometry, the connection between 
the availability of "which-path" information and the loss of coherence has 
been thoroughly examined.E3 How do these results generalize to condensates? 




5. ADDITIONAL DIRECTIONS 
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5.1. Josephson Effects 

In the discussion so far the coupling drive was operated in the very 
short-pulsed mode. If instead it is applied as a weak, continuous drive, it 
has correspondingly smaller bandwidth, such that it is resonant in certain 
regions of the condensate, and off-resonant elsewhere. Population transfer 
can be constrained to occur only within a particular spatial region, for in- 
stance right at the overlap region, and can show oscillatory but nonlinear 
behavior. We are performing experiments!! in this regime, and are explor- 
ing the close analogyS'll between the resulting population transfer and a 
nonlinear Josephson junction. 



5.2. Number States 



In experiments to date, the relative phase appears random if measured 
more than 150 ms after the condensates are split. This could be due to 
intrinsic decohering mechanisms, or it could be due to technical instabilities. 
It is diverting to contemplate, however, creating a condensate pair for which 
from the very beginning the relative phase is intrinsically unpredictable. 
The conjugate variable to the relative phase is the number difference AiV = 
Ni — N2. For experiments with relatively small total N, say less than 2000, 
it may be possible to create states for which N\ ~ N2 ~ N/2, with the 
uncertainty in AiV less than plus or minus 1 atom. Such a state, analogous 
to a highly squeezed state of light, will have a completely undefined relative 
phase and could be useful for performing sub-shotnoise spectroscopy,!! and 
illustrates some instructive paradoxes of quantum measurement theory^ 
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